Superfluidity of spin-1 bosons in optical lattices 

Rong Cheng and J.-Q. Liang 

Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China 
In this paper we show that the superfluidity of cold spin-1 Bose atoms of 

weak interactions in an optical lattice can be realized according to the excita- 
tion energy spectrum which is derived by means of Bogliubov transformation. 
The characteristic of the superfluid-phase spectrum is explained explicitly in 
terms of the nonvanishing critical velocity, i.e., the Landau criterion. It is 
observed that critical velocities of superfluid are different for three spin com- 
ponents and, moreover, can be controlled by adjusting the lattice parameters 
in practical experiments to detect the superfluid phase. 
PACS numbers: 03.75.Lm, 03.75.Mn, 32.80.Pj 

I. INTRODUCTION 

Originally discovered in the system of liquid helium and later in the context of supercon- 
ductors, superfluidity is a hallmark property of interacting quantum fluids and encompasses 
a whole class of fundamental phenomena such as the absence of viscosity, persistent cur- 
rents and quantized vortices. With the achievement of Bose-Einstein condensation (BEC) 
in alkali-metal atoms, the weakly interacting Bose gases have served as an idealized model 
of the superfluid [1] and a test ground of macroscopic quantum effects at low temperatures. 
Recently, rapid advances of experimental techniques in optical traps [2-4] open up a prospect 
for the study of the superfluidity of BEC trapped in periodic potentials, which has attracted 
fast growing interests both experimentally and theoretically [5,6]. The main reason is that 
the optical lattices possess controllable potential depths and lattice constants by adjusting 
the intensity of the laser beams in realistic experiments and moreover the lattice is basically 
defect free. In addition, the great advantage of optical traps is that it liberates the spin 
degree of freedom and provides us an opportunity to test spin-dependent quantum phe- 
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nomena that are absent in the scale-condensate cases. Theoretical studies have predicted 
a variety of novel phenomena of spinor condensates such as fragmented condensation [7], 
skyrmion excitations [8] and propagation of spin waves [8,9]. The quantum phase transition 
from a superfluid to a Mott-insulator (SF-MI) phase in spinor BEC has been observed in 
experiments [4]. In the Mott-insulator phase (MIP), atoms are localized; the particle-number 
fluctuations at each lattice site are suppressed so that there is no phase coherence across the 
lattice. When the tunnel coupling through the interwell barriers becomes large compared 
to the atom-atom interactions, the system undergoes a phase transition into the superfluid 
phase (SFP) in which the atom number per site is random and hence wave function exhibits 
long-range phase coherence. That means one can go from the regime in which the interaction 
energy dominates (high barrier of periodic potential) to the regime where the kinetic energy 
is the leading part (low barrier of periodic potential) by varying the intensity of laser beams 
and vice versa. The SF-MI transition has been also investigated in Refs. [10-14], where the 
Bose-Hubbard model is introduced as the starting point, and the analytic phase-transition 
condition and phase diagram have been obtained by using a perturbation expansion over 
the superfluid order parameter with on-site zero-order energy spectrum which, although 
gives rise to a reasonable description of the MIP, the SFP is not described explicitly. In 
the present paper we study the spinor BEC in a parameter region such that its ground 
state stands deeply in the SFP [14]. The Bogliubov approach is used to obtain the explicit 
excitation energy spectra of weakly interacting spin-1 atoms in an optical lattice and the 
superfluidity is explained explicitly in terms of the energy spectra. The critical velocities 
known as the Landau criterion for the superfluid phase are evaluated for the 23 Na atoms 
and are seen to be realizable in practical experiments. It is also demonstrated that the 
critical velocities are spin component dependent and controllable by adjusting of the lattice 
parameters. Our result may throw light on the experimental observation of the persistent 
atom-current density for the spinor-atom matter waves. 
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II. BOGLIUBOV METHOD AND ENERGY SPECTRUM 



Alkali-metal atoms with nuclear spin / = 3/2, such as 23 Na, 39 K, and 87 Rb, behave at 
low temperatures like simple bosons with a hyperfine spin / = 1. The most general model 
Hamiltonian for the dilute gas of bosonic atoms with hyperfine spin / = 1 trapped in the 
optical potential can be written, in the second- quantization notation, as 

^ = E / d " X ^ PO ("7^ + Vo (X) + V T (A)) & (X) 

+ Y E / d " x ^ (*) $ W (*) & (*) 

a,P,a ,0 

where M is the mass of a single atom; ^ a (X) is the atomic field annihilation op- 
erator associated with atoms in the hyperfine spin state \f = 1,to/ = a) and the 
indices OL,f3,ot,0 label the three spin components [a, (3, a, — — 1, 0, l). Vo (A) = 
K, (sin 2 kXi + sin 2 /cA 2 + sin 2 kX^j is the optical lattice potential formed by laser beams, 
which is assumed to be the same for all three spin components, where k = 2n/\ is the wave 
vector of the laser light with A being the wavelength of the laser light and V is the tunable 
depth of the potential well and, hence, the lattice constant is d — A/2. Vr (A) denotes an 
additional (slowly varying) external trapping potential, e.g., a magnetic trap. The 3x3 
spin matrices F denote the conventional three-dimensional representation (corresponding 
to the spin f = 1) of the angular momentum operator with F%g = {S a ,p-i + $a,p+i) /V%, 
^ap = i ($a,p-i — 8a,p+i) / V%, F*p = a5 a g. The coefficients C and C 2 are related to scat- 
tering lengths ao and a 2 of two colliding bosons with total angular momenta and 2, re- 
spectively, by C = 471-fi 2 (2a 2 + a ) /3M and C 2 = knh 2 (a 2 - a ) /3M. For atoms 23 Na, 
we have a 2 > a , that is C 2 > and the interaction is antiferromagnetic. While, for 87 Rb 
atoms the situation is just opposite that a 2 < ao (this leads to C 2 < 0) and the interaction 
is ferromagnetic [8]. For the periodic potential, the energy eigenstates are Bloch states. We 
can expand the field operator ip a (A) in the Wannier basis, which is a superposition of Bloch 
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states such that 

$ a (X)=Y ( a ai w(X-X i ), (2) 

i 

where w (X — Xj) are the Wannier functions localized in the lattice site i and a ai corresponds 
to the bosonic annihilation operator on the ith lattice site. Also Eq. (2) suggests that atoms 
in different spin states are approximately described by the same coordinate wave function, 
which is seen to be the case when the spin-symmetric interaction is strong compared with the 
asymmetric part, i.e., |C | 3> |C 2 | [11]. This is relevant to experimental conditions for 23 Na 
and 87 Rb atoms. Using Eq. (2), the general Hamiltonian (1) reduces to the Bose-Hubbard 
Hamiltonian 

H = -J YY^^ + Y Y £ Afi™ 

+ \ u ° Y Y alAwai ( 3 ) 

1 a,/3,a',/3' 

Here the first term in Eq. (3) describes the strength of the spin-symmetric tunneling, which is 
characterized by the hopping matrix element between adjacent sites i and j. The tunneling 
constant J = - J d 3 Xw* (X - X { ) [-V 2 /2M + V (X)] w (X - Xj) depends exponentially 
on the depth of potential well Vo and can be varied experimentally by several orders of mag- 
nitude. The second term denotes energy offset of the ith lattice site due to the external con- 
finement of the atoms, where the parameter £i = j d 3 Xw* (X — Xj) Vt {X) w (X — X 4 ) is as- 
sumed to be of the same value e for all lattice sites in the present paper. The third and fourth 
terms characterize the repulsion interaction between two atoms on a single lattice site, which 
is quantified by the on-site interaction matrix element C/ (2) — Co(2) / \ w (X — Xi)\ 4 d 3 X. In 
our case the interaction energy is very well determined by the single parameters Uq and U2, 
due to the short range of the interactions, which is smaller than the lattice spacing. 

Now we use Bogliubov approach to diagonalize the Hamiltonian and obtain the excitation 
energy spectrum of spin-1 bosons with weak interaction [14] in an optical lattice and hence 
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to study the SFP property explicitly. To this end we firstly express the site space operator 
a ai in terms of the wave-vector space operator ak, a as 

aai = ~7K ak ^ e%k ^ 

C-t _ 1 V p-ik-Xi 
a ai ~ v ^vJ Z^k u k,a° ' 

where N s is the total number of the lattice sites and Xj is the coordinate of site i. The 
wave vector k runs only over the first Brillouin zone. In the tight-binding approximation 
(TBA) if we limit our description to simple cubic lattice and substitute Eq. (4) into the 
Hamiltonian (3), we can obtain 

H = Y Y 6 "k,cflk, a + Hint, 
a k 

H in t = 7^- Y Y ( 5 ) 

S a,P k,p,k',p 

+ ^ Y Y 6 k+ P ,k'+p'^l a ^l,/3 F aa' F p^p',^k',a^ 
a,/3,a ,/?' k,p,k' ,p' 

where e (k) = e — Jz cos {kd) with z the number of nearest neighbors of each site. Since the 
number of atoms condensed in the zero-momentum state is much larger than one, we have 

Y = Y "lA) + 1 ~ Y Na0 = ^° > 1 • ( 6 ) 

a a a 

N a0 is the number of condensed atoms of spin-a component in the zero-momentum state 
and A/o is the total number of condensed atoms. So, we can replace the operator a a o and 
a/ a0 with a " c " number y/N a0 . Thus we have 

N a0 = N a - Y®kcfik,a, (7) 

where N a is the total number of atoms of spin-a component and the term of k ^ is exclusive 
in the wave-vector sum. Moreover, in the interacting part of the Hamiltonian when k ^ 
we regard a\ a , ak, a as the small deviation from the operators of vanishing momentum, thus 
all products of four boson operators are approximated as quadratic form, for example, 

^o, a a{ a a , a a 0ia = N 2 aQ ^N 2 a - 2N a a{,cflk,a, (8) 



a,f3,k=£0 



a,/3,k=/=0 



/3,fc^0 



k,j3 a k,P 



N X z] k,cflk, a 



(9) 



where TV = J2 a N a is the total number of atoms. With the approximation Eqs. (8) and (9) 
the total Hamiltonian (5) can be written as 



H = ^N 2 + N(e-zJ) + ^[(N, - N_,) 2 + 2N {y/N~ 1 + ^/N~l) 2 } + X 



2N< 



k^O 



2N. 



a 

J- X V^o V / ^o(a t fc,aa-fc,/3 + ^k, a a- k ,p + 2a\ ia a k ,p) + (2^ (al,i«-fe,-i + a-fc A-i) + X! 



a,(3 7=± 

+a fe)7 a_ fe;7 + 2aj^ j7 a fcj7 ) + 2 ^ v^^V^^(«!, 7 aL fe)0 + a_ kn a kfl + a^a^.o + a^ a fc)7 + 2a£ j0 a fe _ 7 + 2< 



7=±1 



+2a//Vi^ v / ALi (aj[. )0 aL fe)0 - a^ 1 a^ A , _ 1 + a kfi a_ kfi - a fc _ia_ M - aj^aj. _i - a[ _{a k ,i ~ 2a^ a fci0 ) 

Note that the ratio U2/U0 is proportional to the ratio of C2/C0 for all lattice geometries 
and hence U 2 /Uq is small enough in general. Therefore the spin-asymmetric part of the 
interaction is much smaller than the spin-symmetric one. 

The Hamiltonian (10) is quadratic in the operators a k , a , Q>- ka an< ^ can De diagonalized 
by the linear transformation 



a k,a = U k,ab kta — ffc,a^-fe,a, 



(11) 



known as the Bogoliubov transformation. This transformation introduces a new set of 
operators b k>a and b\ a to which we impose the same Bose-operator commutation relations 
bka, b\, = S k k '5 a/ 3 . It is easy to check that the commutation relations are fulfilled if the 



parameters u k a and v k a satisfy the relation 



2 2 
u k,a ~ V k,a = 1) 



(12) 



where the auxiliary parameters u kja and v k>a are to be chosen in order to have the vanishing 
coefficients of the nondiagonal terms b\ a tf_ ka and b k ^ a b_ k ^ a in the Hamiltonian (10) (see the 



Appendix for detail). In virtue of the Bogliubov transformation (11), we finally obtain the 
diagonalized Hamiltonian as 



H = E c + ^ E k , a , a b\ a b k , a (13) 

with 

E c = ^U Nn 2 + N(e-zJ) + ^-[(^i - A_ x ) 2 + 2N Q (y/N' 1 + ^/N^) 2 }, 
where the energy spectra E k ^ a (& = 0, ±1) of the quasiparticle are given by 

£fc, 7 ,7 = \j^k (£k + 2U n 7 + 2U 2 n y ), (14) 



E k ,o,o = \j e k {e k + 2U n + 4U 2y /n^ y/n^) , (15) 

where 7 = ±1 and 

£k(k^o) = zJ[l - cos (kd)] . 

The symbol n a = N a /N s represents the average atom number of spin-a component per 
lattice site. In the experiments of Ref. [4] the number of atoms per lattice site is shown to 
be around 1 — 3. 

III. CRITICAL VELOCITY OF SUPERFLUID 

The energy spectra Eqs. (14) and (15) are typical for the superfluid. To this end we look 
at the dispersion relations of energy spectra E k ^ a for the limit case k — > 

E kaa ~ [zJd 2 (U + U 2 ) n 7 ] 1/2 k, (16) 

E kfi , ~ [zJd 2 (U n + 2E/ 2v /nZT v ^)] 1/2 k. (17) 

The linear wave-vector dependence of the excitation spectra E k ^ a is the characteristic of 
the superfluid which gives rise to critical velocities of superfluid found as 



V 



( 



dE k . 
dk 



.7:7 



) 



h 



1 



[zJd 2 (U + U 2 ) n. 



"i. 



1/2 



(18) 



( 



dk 




1 

h 



[zJd 2 (U n + 2U 2 y/K~[^[)] 



1/2 



(19) 



which reduce to the critical velocity of superfluid given in Ref. [13] for the spin-zero case 
when U 2 = 0, where 1/H is dimension correction. The nonvanishing velocity is nothing 
but the Landau criterion for the superfluid phase. As seen from the above formulas (18) 
and (19), whether there exist critical velocities of superfluid or not depend on appropriate 
values of J and C/q(2) which are related to the Wannier functions determined essentially 
by the potential of optical lattice. Thus, J and C/ (2) can be controlled dependently by 
adjusting the laser parameters. Since the spin-asymmetric interaction U 2 is typically one to 
two orders of magnitude less than the spin-symmetric interaction U , it can ensure that the 
nonvanishing v s , a (a = 0, ±1) exist, whether for the antiferromagnetic interaction (U 2 > 0) 
or for the case of ferromagnetic interaction (U 2 < 0). 

Certainly, it is important to see whether or not the superfluid phase can be realized 
practically with the recent progress of experiments on the confinement of atoms in the light- 
induced trap. To see this we evaluate the values of critical velocities of superfluid adopting 
the typical experimental data in Ref. [15] for a spin-1 condensate of 23 Na atoms in the 
optical lattice created by three perpendicular standing laser beams with A = 985 nm. The 
scattering lengths for 23 Na atoms are ao= (46±5)a B and a 2 = (52 ± 5)a B , where a# is the 
Bohr radius (corresponding to a ratio value U 2 /Uq = 0.04). The valid condition of the 
Bogliubov approach that U « J (for example, U/J « 0.1) can be fulfilled in a region 
of barrier-height values of the optical lattice potential V from to two or three times E R , 
where E R is the recoil energy. It turns out that the magnitude of the critical velocities of 
the superfluid is the order of mm/s which is seen to be in the range of experimental values 



The critical velocities of superfluid v SyOL different for three spin components are functions 
of densities. Therefore, the critical velocities of superfluid can be detected experimentally by 



[15]. 
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counting the atom-number populations. Moreover, the component-dependent velocities may 
imply component separation of spinor BEC in an optical lattice similar to the experimentally 
observed component separation in a binary mixture of BECs [16] since one can control the 
subsequent time evolution of the mixture of a three-condensate system and detect the relative 
motions of the three components which tend to preserve the density profiles, respectively. 
In particular, we may consider a condensate of spin-polarized atoms which are all in the 
state of spin component a = at initial time t — 0, i.e., |-0 (0)) = |0, iV , 0). In this case 
a pair of atoms in the a = state can be excited into the a — ±1 states respectively. 
Thus after a time t c the number of atoms of the a = component becomes steady such 
that N (t c ) = N /2 [11] and the number of atoms of a — ±1 components is about half of 
Nq (t c ), i.e., N_i (t c ) = Ni (t c ) = N (t c ) /2 « N /A. Consequently, the critical velocities of 
superfluid for a = — 1 (f s ,-i) and a — 1 (f Sj i) are equal and different from that for the 
a = component (f s ,o)- There exists a simple relation that f s ,-i = v Sj i = v S fi/V2 for the 
case considered and therefore the atoms of the a = component can be separated from the 
mixture of BECs. 

IV. CONCLUSION 

The energy-band structure of excitation spectra derived in terms of Bogliubov transfor- 
mation for spin-1 cold bosons in an optical lattice is shown to be typical for the superfluid 
phase from the viewpoint of the Landau criterion. Our observation is that the critical veloci- 
ties of the superfluid flow are spin- component dependent and can be controlled by adjusting 
the laser lights that form the optical lattice. The theoretical values of critical velocities ob- 
tained are in agreement with experimental observations and possible experiments to detect 
the superfluid phase are also discussed. 
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VI. APPENDIX 



(Al) 



The inverse transformation of Eq. (11) is 

bk,a = u k,a a k,a "I" v k,oS^— k ,ai 
^k,a = Uk ^k,a + V k , a a-k,a- 

The Hamiltonian Eq. (10) can be written in terms of the quasiboson operators bk, a and b\ a 
as 



H — E c + Hi + H 2 , 



(A2) 



where 



E, 



U 



c ~ 2N S N2 + N ( £ - zJ ) + ^[(^i - + 2N o(\fN~i + jN~i?] 



(A3) 



is a constant. Hi and H 2 denote the diagonal and off-diagonal parts, respectively, with 

fc^O VL\s s J \ s s 



Uo 



10 M fci o"Ufc,0 



- U ° AT J_ U * M \( 2 ^ 2 v /2C/ ^ 2C/ 2 



(u fc)7 t; fc _ 7 + u fc _ 7 u fc)7 ) + - a/A^a/ NZ^(u k ^u k _ 7 + v fe _ 7 
f T / t/o C^ 2 \ / / 2C/ 2 , > 1 . . / U U 2 \ , r 

+ 2^ | [{W + WJ V ^OV^OO + V^-tOV ^00 J («fc )7 «fc,0 + U fc)0 Ufc i7 ) - ^— + — J V^toV j 



,0^,7 

)}®A0 + &lA,7)) 



+ const, 



2iV, 



^0 ^2 / / \ 

Ek + W N ° + 2 A TV iV -10V / ^10 I Mfc,0^fc,0 



(4,0' 
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+ 
+ 



Y u 


u 2 \ 


[\2N S 


N S J 



(Mfc,l M fe -1 + ^Jfc,l^Jfc,-l) 



N s N s 



7 

+ £ 



Uo U, 
2N. 



s 1 v s 

- , u r. 



y/N_ 10 y/N^u kil v ki 



u kn v ka \ (6fc )7 6L fci7 + b ka b_ ka ) 



S f \ S S / s 



u, 



,7&-fc,o) 



C/ Q C/ 2 \ [/ 2 , ^~ 

AT + AT ) V^V^ + 2^ v^V^ 



In order to eliminate the off-diagonal part i?2 we require that the coefficients of all terms 



^l,a^Lfc Q ' an d b k ^ a b_ k a > vanish. In view of condition (12), it is easy to introduce a set of 
parameters (p k>a such that 



U k , a = COsh0 fcjQ , 

v k)a = sinh 4> k>a 



(A6) 



for the convenience of calculation. Conditions (12) and (A6) lead to the useful relations 

2ll k aVk,a 



tanh20 feja = 



2 ' 
k.a ^ u k,a 



cosh (20 M ) = u\ a + v\ a = =, 

y 1 - tanh 2 20 fcjQ 

with which the Hamiltonian (10) can be finally reduced to the diagonal form as given in Eq. 
(13). 
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